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Abstract 

We find classically stable solitons (instantons) in odd (even) dimensional scalar non- 
commutative field theories whose scalar potential, V(4>), has at least two minima. These 
solutions are bubbles of the false vacuum whose size is set by the scale of noncommu- 
tativity. Our construction uses the correspondence between non-commutative fields and 
operators on a single particle Hilbert space. In the case of noncommutative gauge theories 
we note that expanding around a simple solution shifts away the kinetic term and results 
in a purely quartic action with linearly realised gauge symmetries. 
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1. Introduction 

Quantum field theory on a noncommutative space is of interest for a variety of rea- 
sons. It appears to be a self-consistent deformation of the highly constrained structure 
of local quantum field theory. Noncommutative field theories are nonlocal; unraveling the 
consequences of the breakdown of locality at short distances may help understanding non- 
locality in quantum gravity. The discovery of noncommutative quantum field theory in a 
limit of string theory [III] provides new inroads to the subject. 



Perturbative aspects of noncommutative field theories have been analyzed in J2 - 27 ] . 
This study has thrown up some evidence for the renormalizability of a class of noncom- 
mutative field theories, and has revealed an intriguing mixing of the UV and IR jnjj in 
these theories. In this paper we will construct localized classical solutions in some simple 
noncommutative field theories. We expect these objects to play a role in the quantum 
dynamics of the theory. 
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We first consider a scalar field with a polynomial potential. A scaling argument due to 



Derrick |28| shows that, in the commutative case, solitonic solutions do not exist in more 
than 1 + 1 dimensions, as the energy of any field configuration can always be lowered by 
shrinking. Perhaps surprisingly, for sufficiently large noncommutativity parameter 9, we 
will find classically stable solitons in any theory with a scalar potential with more than 
one local minimum. These solitons are asymptotic to the true vacuum, and reach a second 
(possibly false) vacuum in their core. They cannot decay simply by shrinking to zero 
size because sharply peaked field configurations have high energies in noncommutative 
field theories. These solitons are metastable in the quantum theory, but by adjusting 
parameters in the scalar potential, their lifetime can be made arbitrarily long while their 
mass is kept fixed. Solutions are found corresponding to solitons in 21 + 1 dimensions or 
instantons in 21 dimensions for any I. 

Our construction of these solutions exploits the connection between non-commutative 
fields and operators in single particle quantum mechanics. Under this correspondence, the 
* product maps onto usual operator multiplication, and the equation of motion translates 
into algebraic operator equations. The noncommutative scalar action can be rewritten as 
the trace over operators (which can be regarded as oo x oo matrices). This leads to a 
connection between noncommutative field theories, and zero dimensional matrix models. 

Next we consider noncommutative U(N) Yang-Mills theory. When expanded around 
a simple solution of the equations of motion, the action takes the simple quartic form (up 
to constants and topological terms) 

Sym = -r4— / Art^Tr([$ M ,^][$ A ,$ p ]), (1.1) 

where $^ are N x N hermitian matrices and all commutators are constructed from the * 
product. Note that the kinetic term has been shifted away! The usual space-time gauge 
symmetries act linearly as unitary transformations on the fields $ M , and the $ M = 
vacuum leaves even local gauge symmetries unbroken. This construction is similar to that 
of [ |29[ , in which the kinetic term of Witten's string field theory action [[50] is shifted away. 
Indeed, our search for such a construction in noncommutative field theory was motivated 
by the tantalizing analogy, noted in [ 15| , between noncommutative field theories and string 
field theories. The existence of the formulation ( |1 . 1| ) of noncommutative gauge theories 
strengthens the analogy. We also reproduce, as an illustration, the U (1) instanton solutions 
of IBl . 
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Rewriting noncommutative fields as the large N limit of matrices, (LI) is closely 
related to the IKKT matrix theory |32] . Indeed, our construction is essentially equivalent 
to that presented by Aoki et. al. |TT| in this context. Related observations are also made 
in IMSil- 

This paper is organized as follows. In section 2 we describe the action for noncom- 
mutative scalar field theory. In section 3 we consider the limit, 9 — > oo, in which the 
equations simplify considerably. The general solution can be found exactly and is given 
in terms of quantum mechanical projection operators. In Section 4 we show that there 
are stable solitons in this limit, as long as the potential has at least two local minima. In 
section 5 we argue that there are stable solitons at large but finite 9 which can be con- 
structed perturbatively in 9~ x . In section 6 we turn to the noncommutative gauge theory 
where the purely quartic action is constructed. The U(l) instanton solution of [|3]J is also 
reproduced. In an Appendix we give an explicit construction, of the leading i correction 
to the simplest stable soliton of the scalar field theory. 



2. The Noncommutative Scalar Action 

Consider first a noncommutative field theory of a single scalar </> in (2 + 1) dimensions 
with non-commutativity purely in the spatial directions. The spatial R 2 is parametrized 
by complex coordinates z, z. The energy functional 

E=\ I ' d 2 z + V (</>)) , (2.1) 

where d 2 z = dxdy. (We will comment on the generalization to arbitrary dimensions in 
the appropriate places.) Fields in this non-local action are multiplied using the Moyal star 
product, 

(A-kB)(z,z) = e^ a ' a ''- 9 ' ,d ^A{z,z)B{z , ,S')\ z=gt . (2.2) 

Note that in the quadratic part of the action, the star product reduces to the usual product. 

We seek finite energy (localized) solitons of (|2.1| ). These can also be interpreted as 
finite action instantons in the two-dimensional euclidean theory. We will, however, refer 
to the solutions as solitons in the following. 



Since no solutions exist in the commutative limit 9 = |2S[] , we begin our search in the 
limit of large noncommutativity, 9—>oo. It is useful to non-dimensionalize the coordinates 
z^z\^9, z^z\l~9. As a result, the ★ product will henceforth have no 9; i.e. it will be given 
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by (|2.2| ) with 6 = 1. Written in rescaled coordinates, the dependence on 6 in the energy is 
entirely in front of the potential term: 

E = ^J + (2-3) 

In the limit 6—>oo, with V held fixed, the kinetic term in ( |2.3| ) is negligible in comparison to 
V((p), at least for field configurations varying over sizes of order one in our new coordinates. 

Our considerations apply to generic potentials V(4>), but we will, for definiteness, 
mostly discuss those of polynomial form 

y(0) = ImV 2 + E^'- ( 2 - 4 ) 
1 3=3 3 

We have, of course, abbreviated 

(f) J = (f) -k (f) -k ■ ■ ■ -k (f). 



3. Scalar Solitons in the 6 = oo Limit 

After neglecting the kinetic term, the energy 

E = L J d^V^), (3.1) 
is extremised by solving the equation 

dV 

«*=°- (3 - 2) 

For instance, ( |3~2| ) is 

m 2 4> + b 3 (p ★ <j> = (3.3) 

for a cubic potential and 

m 2 (fi + b 3 (j) * (p + b 4 (p * (p * (p = (3.4) 

for a quartic potential. 

If V(<p) were the potential in a commutative scalar field theory, the only solutions to 
( p.2|) would be the constant configurations 



</> = A l5 (3.5) 

where A, G {Ai, A2, ■ ■ ■ , A^} are the various real extrema of the function V(cc)0. As we shall 
see below, the derivatives in the definition of the star product allow for more interesting 



solutions of ( |3.2| ). 



1 For V((f>) as in ( |2.4j ), Ai are the real roots of the equation m 2 x + Y2 r ~3 ^i x ^ 1 = 0- 



j = 3 
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3.1. A Simple Nontrivial Solution 



A non-trivial solution to (|3.2|) can easily be constructed. Given a function <po(x) that 
obeys 

(<Po*<Po) (x) = (/>o(x), (3.6) 

it follows by iteration that <Pq(x) = </>o(x)J and that f (a(po(x)) = f(a)(po(x) (fields in 
/ are multiplied using the star product). In particular, Ai^o(x) solves ( |3.2|) when Aj is 
an extremum of V(x). Thus, in order to find a solution of ( |3.2|) , it is sufficient to find 
a function that squares to itself under the star product. We proceed to construct such a 
function below. 

If we take the ordinary product of a smooth function of width A with itself, the 
spatial size of the function shrinks to a fraction of A, which is why non-constant functions 
never square to themselves! The non-locality of the star product, however, introduces an 
additional effect, adding roughlyi i to the width of the product. This makes it possible 
for a lump of approximately unit size to square to itself under the star product. 

Consider a gaussian packet of the form 

1 _ - 2 
ipA{r) = —r^e a*, 

with radial width A (here r 2 = x 2 + y 2 ). The star product of ijj A with itself is easily 
computed by passing to momentum space, 

Mk) = I e lk - x i(; A (x)d 2 x = e"^, (3.7) 



(Va * ^a) (p) = J d 2 kip A (k)ip A (p - fc)e*' 

2ttA 2 

Therefore 



A 2 + ^ 



(3.f 



(3.9) 



In particular^, when A 2 = 1, the gaussian squares to itself (up to a factor of 2%). That is, 

<f>o{x) = 2irip 1 (x) = 2e~ r2 (3.10) 
solves ( p.6|) and Xi4>o(x) solves ( |3.2|) . 



This equation and its solution has also appeared in earlier work involving the Moyal Product. 
See H||39]. 

3 The added width is actually ~ K, the typical momentum in the Fourier transform of the 
function. For a function of size A with no oscillations, K ~ ^. For a function of size A with n 
oscillations, K « ^. 



We note in passing that in the limit A^O, (|3.9|) reduces to 5 2 (x) * S 2 (x) 
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3.2. The General Solution 

In order to find all solutions of ( |3.2|) we will exploit the connection between Moyal 
products and quantization. Given a C°° function f(q,p) on R 2 (thought of as the phase 
space of a one-dimensional particle) , there is a prescription which uniquely assigns to it an 
operator Of(q,p), acting on the corresponding single particle quantum mechanical Hilbert 
space, TC. It is convenient for our purposes to choose the Weyl or symmetric ordering 
prescription 

Of(q,p) = J^yJ d 2 kf(k)e-< k ^ +k ^ 1 (3.11) 



where 



and 



f(k) = J d 2 xe l ^ q+k o p) f(q,p), (3.12) 



[q,p\=i. (3.13) 
With this prescription, it may be verified that 

i- J dpdqf(q,p)=Tr n O f , (3.14) 

and that the Moyal product of functions is isomorphic to ordinary operator multiplication 

O f -O g = O f + g . (3.15) 

In order to solve any algebraic equation involving the star product, it is thus sufficient 
to determine all operator solutions to the equation in TC. The functions on phase space 
corresponding to each of these operators may then be read off from (|3.11| ). We will now 



employ this procedure to find all solutions of (|3.2|) . 

As noted above, any solution to (|3.6|) may be rescaled into a solution of (|3.2|) . Par- 
ticular solutions of ( |3T2] ) may thus be obtained by constructing operators in TC that obey 
(|3~6|), i.e. O 2 ^ = O^. This equation is solved by any projection operator in TC. TC possesses 
an infinite number of projection operators, which can be classified by the dimension of the 
subspace they project onto. Each class contains a large continuous infinity of operators, 
each of which, upon rescaling, yields a solution to (|3~2|). 



The most general solution to ( |3.2| ) hence takes the form 

Q = ^a 3 Pj (3.16) 



3 
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where {Pj} are mutually orthogonal projection operators onto one dimensional subspaces, 
with cij taking values in the set {Aj} of real extrema of V(x). 

In order to obtain the functions in space corresponding to the solutions ( |3.16|) , it is 



convenient to choose a particular basis in TC. Let \n) represent the energy eigenstates of the 
one dimensional harmonic oscillator whose creation and annihilation operators are defined 
by 

a= ^T ~JT' (3 ' 17) 

Note that a\n) = \fn\n — 1) and a^\n) — vn + l\n + 1). Any operator may be written as 
a linear combination of the basis operators |m)(n|'s, which, in turn, may be expressed in 
terms of a and as 

fm n 

\m)(n\ =: e~ a a -= : (3.18) 
y/ml Vn! 

where double dots denote normal ordering. 

We will first describe operators of the form (|3.16| ) that correspond to radially sym- 

metric functions in space. As a T a ~ operators corresponding to radially symmetric 

wavefunctions are functions of a^a. From (|3.18| ), the only such operators are linear combi- 



nations of the diagonal projection operators \n)(n\ = A : a^ n e~ a ^ a a n :. Hence all radially 

symmetric solutions of (|3.2| ) correspond to operators of the form O = a n \n){n\, where 
the numbers a n can take any values in the set {A^}. 

We now translate these operator solutions back to field space. From the Baker- 
Campbell-Hausdorff formula 



where 

k x + ik y _ k x - ik y 2 _ 

tv z — — y= , l\i Z — — y= , ft — £tv z tv z . 

Any operator O expressed as a normal ordered function of a and a\ /7v(a,a^), can be 
rewritten in Weyl ordered form as follows. By definition, 

O =: Ma, at) := JL J d 2 kf N {k) : e ~K^+^) . (3 2Q) 

Using ( 3.19 ), (|3.20|) may be rewritten as 

O = — ^ / d 2 kf N (k)e^e-< k " a+k ^). (3.21) 

J 
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Thus, the momentum space function / associated with the operator O, according to the 
rule fl-rnp is 

f(k)=e^f N (k). (3.22) 

For the operator O n = \n)(n\ we find, using ( |3.18| ) and ( |3.2U| ), that the corresponding 
normal ordered function <p^\k) = 2ne^~ L n (^-). ( |3.22| ) then becomes 

\n)(n\ = / d 2 ke^L n (^-)e-< k * a+k ^) (3.23) 

(27TJ y 2 

where L n (x) is the Laguerre polynomial. The field 4> n (x,y) that corresponds to the 
operator O n = \n)(n\ is, therefore, 

<p n (r 2 = x 2 + y 2 ) = -L [ d 2 ke^L n (^-)e-^ = 2(-l) n e - r2 L n (2r 2 ). (3.24) 
{2tt) J 2 

Note that 4>o(r 2 ) is precisely the gaussian solution found in Sec. 3.1. 

In summary, (|3.2| ) has an infinite number of real radial solutions, given by 

oo 

J2 a nMr 2 ) (3.25) 

n=0 

where (f> n (r 2 ) is given by ( |3.24| ) and each a n takes values in {A^}. 

In order to generate all non radially symmetric solutions to ( |3.2| ), we rewrite ( |3.1[ ) in 
operator language, using (|3.14| ) as 

E = ^W(C^). (3.26) 

g 2 

( 3. 26]) is manifestly invariant under unitary transformations of and so has a U(oo) 
global symmetry. In other words, if O is a solution to the equation of motion, so is UOW , 
where U is any unitary operator acting on TC. A general Hermitian operator (one that 
corresponds to a real field 0) may be obtained by acting on a diagonal operator (i.e. an 
operator that corresponds to a radially symmetric field configuration) by an element of the 
U(oo) symmetry group (since any hermitian operator is unitarily diagonalizable) . Thus 
every solution to ( |3.2| ) may be obtained from a radially symmetric solution by means of 
U(oo) symmetry transformations. 

Therefore solutions to ( |3.2| ) consist of disjoint infinite dimensional manifolds labelled 
by the set of eigenvalues of the corresponding operator. Points on the same manifold can 
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be mapped into each other by U(oo) transformations. Each manifold includes several! 
diagonal operators (radially symmetric solutions). We will have more to say about the 
moduli space of these solutions in the next section. 

As all solutions are related to radially symmetric solutions by a symmetry transfor- 
mation, we will mostly discuss only radially symmetric solutions. 

3.3. UV/IR Mixing 

</>o(r 2 ), the Gaussian solution worked out in subsec 3.1, is a lump of unit size centred 
at the origin, as shown in Fig. 1. 




-2-1 12 



Fig. 1: A plot of 4>o( r ) versus r. The solution is a blob centred at the origin. 



4>n(r), at large n, looks quite different (see Fig. 2.). It is a solution of size « y/n 
that undergoes n oscillationsi in that interval, with oscillation period oc -7=. </> n (r 2 ) thus 
receives significant contributions from momenta up to yfn in momentum space. These 
solutions exemplify the UV-IR mixing pointed out in oscillations with frequency yfn 
produce an object of size y/n (instead of -^=) in a noncommutative theory. 



5 Distinct diagonal operators having the same eigenvalues lie on the same manifold, being 
related by the "Weyl" subgroup of U (00) that permutes eigenvalues. 

6 Using asymptotic formulae for Laguerre polynomials we find 



4>n{r) 



2{-l) n 

2(-l)« 

— - — —r cos( 

(277 2 r 2 ) 4 
2(-2r 2 )" -r S 



r < 



2n2r - 2.) 



< r < V2n) 



r » \/2n 
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Fig. 2: A plot of <p3o(r) versus r. 
3-4- Generalization to Higher Dimensions 

All considerations of the preceding subsections may easily be generalized to higher 
dimensions. Consider a scalar field theory in 21 + 1 dimensions with non-commutativity 
only in the spatial directions. By a choice of axes, the 21 x 21 dimensional noncommutativity 
matrix O may always be brought into block diagonal form. In other words, it is possible 
to choose spatial coordinates Z{,Zj (i,j = in terms of which the non-commutativity 

matrix = ^^j, Qij = G^j = 0., As before we consider the limit where 9i are uniformly 
taken to oo and non-dimensionalize Zi^Zi\f6l. As in the previous subsections, the kinetic 
term in the action may be dropped in this limit. Solutions to the equations of motion 
( ^2|) are once again in correspondence with operator solutions to the same equations; the 
operators in question now acting on H x 7i x ... x 7i, I copies of the Hilbert space of the 
previous subsection. The general solution to (|3.2|) once again takes the form ( |3.16|) in terms 



of projection operators on this space. As in the previous subsection, the general solution 
may be obtained from diagonal solutions via U(oo) rotations. Diagonal solutions to ( |3.2| ) 
are given by 

O = ^2a,fi\n){n\ <-> ^cifi JJ^n^N 2 ), (3.27) 

n n i 

where n is shorthand for the set of quantum numbers {ni} for the I dimensional oscillator 
and <j) ni are defined in ( |3.24j) . As in ( |3.25| ), the coefficients take values in {A^}. A subset 



of the solutions ( |3.27| ) are actually invariant under SO (21) rotations and can be written in 
terms of associated Laguerre polynomials. These are displayed in Sec. A. 3 of the Appendix. 

In summary, in the limit of maximal noncommutativity, the construction of solitons 
in two spatial dimensions generalizes almost trivially to every even spatial dimension. 
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4. Stability and Moduli Space at 9 = oo 

In this section we study the stability of the solitons constructed in the previous section. 
We will also describe the moduli space of stable solitons. 

4-1. Stability at 9 = 00 

We wish to examine the stability of the radial solution 



(4.1) 



n=0 



to small fluctuations. Since any U(oo) rotation does not change the energy of our solution 
( [4~1~1) , it is sufficient to study the stability of (|4.1|) to radially symmetric fluctuations. These 
are most conveniently parameterized as deformations of the eigenvalues. The energy for 
an arbitrary radially symmetric state <p(r 2 ) = 5^^L c n <p n (r 2 ) is 



E 
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n=0 



The solution c n = X an is manifestly an extremum of S, as, by definition, X ai are extrema 
of the function V(x). Clearly ( fl.l|) is a local minimum of the energy (and so a stable 
solution) if, and only if, A a „ is a local minimum of V(x) for all < n < 00. 

As an example consider the cubic potential of Fig. 3. with a maximum at A = —1. In 
this case, all X an in ( |4.1|) are either zero or -1. The only stable solution is that for which all 
X an = 0, i.e. the vacuum. The solution — </> (r 2 ), for instance, is unstable, as the energy of 
this field configuration is decreased by scaling this solution by a constant near unity. This 
instability shows up as a negative eigenvalue of the quadratic form for fluctuations about 
— 4>o( x )', the corresponding eigenmode 5<fio is oc 4>q. 




Fig. 3: The <p 3 theory with an unstable extremum 
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On the other hand the field theory with V(4>) (say, for a quartic potential) graphed 
in Fig. 4 has stable solitons; these are solutions of the form 4>(r 2 ) = Y^=o ^c n 4>n(r 2 ) with 
A c „ taking the values of the minima - or A « —1.4 for all n. In particular \(fio(r 2 ) is a 
stable solution, manifestly stable to rescalings. Again, one may check that the quadratic 
form for fluctuations about \4>o(r 2 ) is positive. In particular, 6(f) oc 0o is an eigenmode of 
this quadratic form with positive eigenvalue. 
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Fig. 4: A (p 4 potential with two minima. 

The stability of <f)(r 2 ) = \(po(r 2 ) in the previous example may qualitatively be un- 
derstood as follows. <pQ is a Gaussian of height 2 A. Far away from the origin, <po{x) = 0, 
but near x = 0, 4>o(x) is in the vicinity of the second vacuum. In other words, the static 
solution corresponds to a bubble of the "false" vacuum. The area of the bubble is of order 
one (or 9 in our original coordinates), the non-commutativity scale. In a commutative 
theory such a bubble would decay by shrinking to zero size. Noncommutativity prevents 
the bubble from shrinking to a spatial size smaller than \f9. In order to decay, 4>q actually 
has to scale to zero - but that process involves going over the hump in the potential and 
so is classically forbidden. 



p. 8 
/ 0.6 




/ 0.4 

/ 0.2 




-2 -1 


1 2 



Fig. 5: Profile of the Gaussian soliton with a false vacuum region (above 
the horizontal bar) of radius 1. 

The energy of this soliton is proportional to the vacuum energy density at the 
'false' vacuum times the volume of the soliton 9. It is remarkable that the energy of the 
soliton is completely insensitive to the value of the scalar potential at any point except 
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(j) = A. Thus the mass of the soliton is unchanged if the height of the barrier in V((p) 
(between </> = X and </> = 0, see Fig. 4.) is taken to infinity while V(X) is kept fixed. 
This is true even though 4>o(r), the solitonic field configuration corresponding to A|0)(0|, 
decreases continuously from </> = 2Aatr = 0to</> = 0atr = oo! 

Consider a 2+1 dimensional scalar theory, noncommutative only in space, at infinite 
9. Using the correspondence between functions and operators (matrices) described in the 
previous section, the noncommutative scalar field theory is equivalent to the matrix quan- 
tum mechanics of an N x N hermitian matrix H, at infinite N, with the usual relativistic 
kinetic term Tr (dtH) , and a potential Tr (V(H)). The amplitude for an eigenvalue of H 
to tunnel from A to is exponentially suppressed by the area under the potential barrier 
in Fig. 4., and goes to zero as this area is taken to infinity. Thus the finite mass soliton 
A 1 0} (0 1 is stable, even quantum mechanically, in this limit. 

The U(oo) symmetry of ( |3.1j ) is spontaneously broken by every nonzero solution, (p(x), 
of ( |3.2|) . As a consequence, every solution has a number of exact zero modes (Goldstone 
modes) corresponding to small displacements about <p(x) on the manifold of solutions. As 
Rnm — \n){m\ + \m)(n\ and S nm = i(\n)(m\ — \m)(n\) are the generators of U(oo), these 
zero modes are given by the nonzero elements of S(f> oc [R nm , <f>], [S nm , </>]. 

The U(oo) group of symmetry transformations that generates these zero modes is 
certainly not manifest (at least to the untrained eye) in the energy written in coordinate 
space in the form ( |3.1p . In addition to the two translations, ( |3.1| ) possesses three manifest 
local symmetries, corresponding to a linear change of the coordinates x,y by an SL(2,R) 
matrix. The remaining U(oo) transformations act non-locally on <p(x,y), according to 
(f)'(x,y) = (U -ktfi-kW) (x,y) where U(x,y) is any function that obeys U ★ U* = 1 (such 
functions correspond to U(oo) operators under the map ( J3.11| )). 

All arguments in this subsection may be applied (after straightforward generalizations) 
to higher dimensional solitons. 

4-2. Multi Solitons 

In this subsection we will qualitatively describe a part of the moduli space of stable 
solitons (at 9 = oo) in the simple case of the potential graphed in Fig. 4 with a single 
non-zero minimum at (j) = A. 

The stable solitons can be characterized by their 'level' (number of A eigenvalues) . All 
stable level one solitons correspond to operators of the form 

XU\0)(0\U^ (4.2) 
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where U is a unitary operator. As mentioned above, the set of level one solitons span an 
infinite dimensional manifold parameterized by U(N)/U(N — 1) (for N = oo). 

The soliton looks very different at different points on the manifold. U = I in ( |4.2| ) 
corresponds to the gaussian blob of Fig. 1. If U happens to be a unitary transformation 
that maps |0) to \m), for large m, the corresponding wave function is qualitatively similar 
to that in Fig. 2. When U = e atz ~ az is the generator of translations, the operator in Q4.2|) , 

A|z)(z|, is proportional to the projection operator onto a gaussian centred around z = 

i i z i 2 t 

-j=(x + iy). (Here \z) = e ~e a z \0) is the usual coherent state.) Again, if U corresponds 
to one of the SL(2, R) operators, we obtain squeezed states; gaussians elongated in the y 
direction and shrunk in the x direction. And so on. 

Turn now to solitons at arbitrary level n. All such solitons may be obtained by acting 

on 

\(4> + 01 H h n -l) 

by arbitrary unitary transformations. The manifold of solutions thus generated is param- 
eterized by j7^v3~) (and has dimension d n ~ 2nN) where N—>oo. 

Notice that d n ~ nd\. This fact has a nice explanation; in a particular limit the 
manifold of level n solutions reduces to n copies of level 1 solitons very far from each 
other. This conclusion follows from the observation that the operator that represents n 
widely separated level one solitons (with centres Zj), for instance 

M = XJ2\zj)(zj\ (4.3) 

j 

is approximately a level n soliton (and exponentially close to a true level n soliton) when 
\zi — Zj\—>-oo for all We demonstrate this explicitly below for the case n = 2. 
Using (z\ — z) = e~ 2 \ z \ , it is easy to check that the kets 

\ Z± ) = !*>*!-*> (4 .4) 

/ 2(1±e -2| 2 p ) ^ > 



are orthogonal. From (|3.16|) we conclude that the projector 



\z)(z\ + | - z)(- z\+e- 2 \ z \ 2 (\z)(-z\ + | - z)(z\) 
(1-e- 



O z = A l + \z_)(z.\) = A " 1 ~ - » p - - 1 - 'j™ (4.5) 



corresponds to a level 2 solution. Up to corrections of order e 2 ' z l , O z is equal to \z)(z\ + 
| —z)(—z\, the superposition of field configurations corresponding to two widely separated 
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level one solitons @. We conclude that a part of the level n moduli space describes n widely 
separated level one solitons. 

We have, so far, worked in the strict limit 9 = oo. The picture developed in this limit 
is qualitatively modified at large but finite 6>, as we will describe in the next section. 



5. Scalar Solitons at Large but Finite 9. 

We have argued that, under certain conditions on V((p), ( |3.2| ) has an infinite number of 
stable solutions. Each solution has an infinite number of exact zero modes, the Goldstone 
modes of the spontaneously broken U(oo) symmetry of ( |3.1| ). 

At finite 9, the kinetic term in (|2.3| ) explicitly breaks this U(oo) symmetry down to 
the Euclidean group in 2 dimensions. Finite 9 effects may thus be expected to 

1. Lift the 9 = oo manifold of solutions to a discrete set of solutions. 

2. Give (positive or negative) masses to the [/(oo) Goldstone bosons about these discrete 
solutions. 

In Appendix A we will argue that, at large enough 9, corresponding to every radially 
symmetric solution s of (|3?2|), there is a radially symmetric saddle point of (|2.3| ), that 
reduces to s as 9^oo. It is likely that these are the only saddle points of Q2.3p. 

Not all these radially symmetric solutions are stable, however. In fact, it might seem 
likely that some of the infinite number of zero modes, at 9 = oo, about each solution s, 



might become tachyonic at finite 9. If this were true, (|2.3| ) would have no classically stable 
extremum at any finite 9, no matter how large. 

We will find that is not the case. In subsection 5.1 below we will argue that any small 
perturbation of ( |3.1| ) must preserve the existence of at least one classically stable level 
one soliton. In subsection 5.2 we will identify this soliton to be the one near the gaussian 

\Mr 2 )- 

5.1. Existence of A Stable Soliton 

For definiteness, through the rest of this section we assume that the potential V(<p) has 
the shape shown in Fig. 4. In particular, it is positive definite. Let the stable extremum 
of V occur at 4> = A and the unstable extremum at 4> = (3, (A < (3 < 0). 



7 It is curious that the kinetic energy of this field configuration is independent of z indicating 
that there is no force between the two solitons even to next leading order in ^. 
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Consider, first, (|3.1| ) i.e. the energy functional in the limit where we neglect the kinetic 
term. We will show that any path in field space leading from the soliton A0o(r 2 ) to the the 
vacuum passes through a point whose energy is larger than ^r-V(fl). Since the energy of 
the stable soliton is ^V(X) < *%rV(P), every path from the soliton to the vacuum must 
pass over a barrier of height 0{\). 

The energy evaluated on an operator A is 

E= 2 4k ( V(A) ) = 2 4±V(c„) (5.1) 

a a n=l 

where c n are the eigenvalues of A. Since V is positive definite, 

E>™V(b), (5.2) 

g 2 

where b is the smallest eigenvalue of A. 

Consider a path in field, or operator space, leading from A0o to the vacuum. At the 
beginning of this path b = A. At its end 6 = 0. Since A < (3 < 0, any smooth path with 
these endpoints must have a point at which b = f3. At that point E > ^rV(P), as was to 
be shown. 

Now include the the kinetic term in ( |3.1| ). Barring singular behaviour, this changes the 
energies of all field configurations by terms of O(^). For large enough the arguments 
of the previous paragraph imply that the field configuration that describe the level one 
soliton at 9 = oo cannot decay to the vacuum. Hence there must exist at least one stable 
soliton near one of the unperturbed level one solutions. In fact, as we will show in the 
next subsection, there is a stable soliton near the gaussian X(po(r 2 ). In the Appendix we 
will present an approximate construction of this solution at large but finite 9. A similar 
argument demonstrates the existence of at least one stable solution at level n. 

5.2. Approximate Description of the Stable Soliton 

All level one solutions to ( |3.2|) take the form XU\0) (0\U* where U is a unitary operator. 
We wish to determine the contribution of the kinetic term to the energy of such an operator. 
The kinetic term in ( p.l|) for an operator A is 

K = ^Tr[a,A][A,a*}. (5.3) 
9 
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Setting A = XU\0)(0\W we find 

= 1 + E 2 ^>o! 2 - 2| £ v/fc+Tt/^o^+Lol 2 - (5-4) 



We expand (|5.4|) to quadratic order in deviations from U = I. Choose U = t/^o for 
i > 1 as the coordinates for this expansion (|C/oo| is determined in terms of Ui as U is 
unitary). To quadratic order in Ui 

?^l = l + 2±k\U k f. (5.5) 

k=2 



As Ui and Ui do not appear in (|5.5|), they parameterize flat directions of K(U) (to 
quadratic order). This was to be expected. Any localized extremum of (|2.1| ) has two exact 
translational zero modes. Infinitesimally, Uqi and its complex conjugates act as derivatives 
on <po(r 2 ), generating these zero modes. Modulo these zero modes, the fluctuation matrix 
about U = 1 is positive definite. 

While K{U) has several critical points other than U = I, it has no further local 
minima. For example, U = U( m \ the unitary transformation that rotates |0)(0| to |m)(m|, 
is an unstable critical point of K(U) for all m. In fact U = U^ is unstable to decay into 
U = I. This may be demonstrated by considering the path in field space where 
\a>) = cosa|0) + sina|m). (|5.3| ) evaluated on such a path is equal to 1 + 2msin 2 a (for 
m > 1; 1 + 2 sin 4 a for m = 1) indicating that the state \m)(m\ can decay to |0)(0|. 

We will now argue that, at large enough 6, the finite 9 saddle point cf>(x, y) of Q2.3Q 
that reduces to A|0)(0| as 9^oo is classically stable. 

Consider the mass matrix for fluctuations about <p(x,y). Since any operator may be 
written as UDU^ where D is diagonal and U unitary, small fluctuations may be decomposed 
into radial (fluctuations of D) and angular ones (fluctuations of U). The mass matrix for 
purely radial fluctuations is O{0) to leading order, and has been shown to be positive 
definite in sec 3.4. The mass matrix for purely angular fluctuations is 0(1) to leading 
order, and has been shown to be positive definite, modulo the two zero modes. Since 
angular modes completely disappear from the potential, mixing between radial and angular 
fluctuations occurs only through the kinetic term, and are also 0(1). These cross terms 
result in corrections to the eigenvalues of the mass matrix only at O(h). Hence, to leading 
order in I, the mass matrix is positive. The two zero modes of the angular mass matrix 
cannot be driven negative by -s corrections as they are exact. 
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A similar argument demonstrates the instability of all other radially symmetric level 
one solitons (those that reduce to A|n)(n| at 6 = oo) at large enough 9. 

The considerations of this subsection may easily be generalized to solitons in 21 spatial 
dimensions, using the higher dimensional analogue of ([5.4]) : 



= 1 + 2 £ kjW ^ - 2£ I E V^T1U^ +J f (5.6) 

j,k j k 



and (Km 



9 2 K(U) 
(2tv) 1 X 2 



We use the notation of ( |3.27| ); k is an / dimensional vector, j runs from 1 to I and i is the 
basis unit vector in the i th direction; in components i n = 8^ n . Notice that K{U) in (|5.7| ) 
is independent of U^ for all i, a consequence of the exact translational invariance in all 21 
spatial directions. 



6. Noncommutative Yang-Mills 

6.1. Quartic Action for the U(l) Theory in Two Dimensions 
Consider the action 

S=j^r- I d 2 z[$,m,*l (6-1) 
where $ is a complex field and 

= (6.2) 

The equation of motion following from ( |6.1| ) is 

[$,[$,$]] =0. (6.3) 

$ can also be viewed as a quantum mechanical operator and $ as it's hermitian conjugate. 
The commutators in ( |6.1| )-( |673[ ) are then ordinary operator commutators, and the integral 
is the trace over the Hilbert space. In the operator representation a simple solution of the 
equation of motion ( |6.3|) is 

$ = a, $ = a f . (6.4) 
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Let us expand around this solution by defining 

$ = a + iA z , $ = a t -iA z . (6.5) 

One then finds, translating back to functions (with \[2z = q + ip, [a, ] = d z and \a) , ] = 
— d z , that 

[$, $] = 1 + id z A- z - id- z A z - [A z , A,} = 1 + %F Z - Z . (6.6) 

The operator representation of ( |6.1[ ) has the manifest U(N = oo) symmetry under which 
= U^U just as in the scalar field theory. Infinitesimally, 

= A], (6.7) 

where = expiA. When gauged, this is just the usual U{1) gauge symmetry of the 
non-commutative theory, 

5A = dA + i[A,A]. 
The equation of motion (|6.3| ) is 

D- Z F Z - Z = 0. (6.8) 

The action (|6.1D is then 

S=~A- I d 2 z(F z - z -t)\ (6.9) 

the standard two dimensional non-commutative U(l) Yang-Mills action up to constants 
and topological terms. 

6.2. The U(N) Theory in 21 Dimensions 
( |BTT|) can be generalized to 

S= -4— / rf 2Z x^ A ^Tr([^,^][$ A ,$ p ]), (6.10) 

where fj,, v = 1, , , 21 and $ M are real N x N matrices. Though we have restricted ourselves 
to a flat euclidean metric, one can generalise the argument below to the Minkowski metric 
as well. 

The equation of motion is 
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(6.11) 



We choose complex coordinates such that Q a i = iS ab -, with a,b = ( |6.11|) has the 

solution 

$ & = a b , $- b = a|, (6.12) 

where [a b , a|] = 5&c- Expanding around this solution with 

$ b = a b + iA b (6.13) 

one finds 



S 



i- / >*(F a5 -e^) 2 . (6.14) 



As before the manifest U(oo)<g>U(N) symmetry corresponds to the non-commutative U(N) 
gauge symmetry. 

6.3. The U(l) Instanton 

The four dimensional non-commutative gauge theory has instanton solutions which 
are deformed versions of the usual non-abelian instantons. In particular, the U(l) non- 
commutative theory also has non-singular finite action saddle points ||31|1 - We exhibit the 
operators $ a corresponding to the simplest such U(l) instanton. 

The operators $ a corresponding to an anti self dual field strength S ab F ab = (a, b = 

1.2) , obey 

[$ 6 ,$c] = 0, S a ~ b [$a,$- b } = 2- (6.15) 

In four dimensions, the operators $ a (fl = 1,2) live in a Hilbert space generated by 
the creation and annihilation operators of a two-dimensional harmonic oscillator (See Sec. 

3.3) . Rather than work in the conventional number basis \ni, 712), it is convenient to work 
in Schwinger's angular momentum basis, 

\j, m) = — =|0, 0), (6.16) 

V(j +m)\ a/(j - m)\ 

with < j < 00, |m| < j. The operators 

J_i_ = a{<22, «/- = a^ai, J z = -{a\a\ — a\a,2) (6-17) 
obey the usual angular momentum algebra. 
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We will find a solution to ( 6.15 ) of the form 



$b = a b 2^(1 + cj)\j, m) (j, m\ = a b + a b ^Cj \ j, m) (j, m\, 

3,m j,m (6.18) 

$b = 4' 

and put it into Hermitian form via a complexified gauge transformation W. 

The ansatz ( |6.18| ) satisfies the holomorphic part of ( ]6.15| ) for any Cj. For a real Cj, 



the only condition comes from the equation Fu = — i*22- Using 



ai,2b',m) = y/j ±m\j - -,m^ -), 
yields the equation jcj = (j + l)c + i. Which has the solution 



1 1, 



(6.19) 



2 



C,- = r, (j > 0). (6.20) 

i(2i + i)' u ; 1 ; 



The complexified gauge transformation 



W = W t = J2 x l-^\j,m) (6.21) 



puts the solution (|6.18|) into Hermitian form for c = — 1. The field strength then takes the 
compact form 

[$ 5 , $ c ] = -5- bc - zF Sc = -5 5c - (J- a) 5c ^ .,. + 1 w 2 . + 1 J i» ™> 0', H- ( 6 - 22 ) 

Here J are the angular momentum generators defined in (|6.171 ) and a, the usual Pauli 
matrices. This solution is exactly the same as the simplest charge one U(l) instanton in 
| 31j| . It may be checked that ^TrF^ = 1. 
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Appendix A. Solutions at Finite 9 

In this appendix we will examine radially symmetric saddle points of ( [2.3| ) at finite 
9. In subsection A.l we study the equation of motion resulting from ( |2.3|) at finite 9, 
and examine the existence of radially symmetric solutions to these equations. In A. 2 we 
concentrate on a particular solution; the one that reduces to the stable soliton A|0)(0| as 
9 is taken to infinity. We present an approximate construction of this soliton at large 
9. In A.3 we briefly comment on the generalization of these results to solitons in higher 
dimensions. 

A.l. The Perturbation Expansion and a Recursion Relation 

The full equation of motion derived from (|2j]) may be written in momentum space as 

3=3 

While the LHS of ( |A.1| ) is independent of 9, the RHS is of order i, and so is a small 
parameter at large 9. For notational convenience, we set -4 = dj and = £■ 
Let 

oo 

Y,cJ n (k 2 ) (A.2) 

n=0 

be a solution to ( |A.1|) . Substituting ([A.2Q into (|A.1|) , using the recurrence relation for 



Laguerre polynomials, and equating coefficients of 4> n (k 2 ), we arrive at the difference equa- 
tions 

r 

c n + S ^2d j c'~ 1 = 2e[nc n -i - (2n + l)c n + (n + l)c n+i ]. (A.3) 

3=3 

We are interested in finite energy solutions to (|3.2| ), i.e. solutions to (|A.3| ) for which 



^l/(c n )<oo. (A.4) 
Since V(0) = 0, (|A.4|) will be satisfied if the c n s approach zero sufficiently fast as n 



approaches infinity. For such a solution, all nonlinear terms in ( A.3 ) may be neglected at 
large enough n. At sufficiently large n, n may also be replaced by a continuous variable u, 
and ( |A.3|) turns into the second order differential equation 

. . d 2 c(u) , . „, 

c(u) = 2eu - V • ( A -5) 
du z 
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( |A.5| ) is the Schroedinger equation for a zero energy state of a particle in a — potential. 
yfe plays the role of Planck's constant, and at small e ( |A.5| ) is easily solved in the WKB 
approximation, yielding 



c(u)=A_u^e v—+A + u*e + v~ (A.6) 

where A± are arbitrary constants. In order that c n tend to zero at large n, A + = 0. Thus, 
for large! n, 

c n ttAn%e~v?. (A. 7) 



( |A.7| ) has an undetermined parameter A, the scale of the solution at large n. As ( |A.3| ) 



is a nonlinear equation, A is not an arbitrary parameter, but is determined to be one of 
a discrete set of values. Given c p and Cp+i, the (p + 1) equations ( |A.3| ) with n = 0- • -p 



overdetermine the p unknowns c n for n < p. The extra equation constrains the scale A, as 
we'll see in the next subsection. 

A. 2. The Gaussian Soliton Corrected 

In this section we present an approximate construction of the stable soliton that 
reduces to the gaussian at infinite 9. Our construction approximates the true solution to 
arbitrary accuracy at small enough e. 

We wish to find a solution of ( |A.3|) such that 

lim c = A (A.8) 



limc m = (A. 9) 



and 



uniformly in m, for m > 1. ( |A.9| ) ensures that, on such a solution, ( |A.3| ) for n > 1 reduces 
to 

c n = 2e[nc n - 1 - (2n + l)c n + (n + l)c n+1 ] (A. 10) 

for small enough e. It is easy to find an explicit solution to ( [A.lOp that obeys ( |A.8| ), ( |A.9| ). 
Consider a function (f>(x, y) that obeys the differential equation 

(-ed 2 + 1)0 = b<p . (A.ll) 



(A.6) is a good approximation when \c n \ <C 1 (so that dropping nonlinear terms in (|A.3|) is 



justified) and c " c c " 1 <C 1, i.e. ne 3> 1 (so that the transition from ( |A.3[ ) to ( |A.5 ) is justified). 
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Expanding <f> in the form 

= J] Cn n (A.12) 

n=0 

and imitating the manipulations of section 4.3, we find that c n s obey ( |A.10|) for n > 1, 
but obey 

c = 2e[ci -c ] + b (A.13) 



instead of ( A. 3 ) (with n = 0). This relation will fix the free parameter b. 



( |A.11|) is easily solved in momentum space 



?W = ^- (A.14) 
Using the explicit forms for 4> n (k) and orthogonality of the Laguerre polynomials we find 



oo 



r„ = h I dx e ~ XLn o {x) . (A.15) 
1 + 2ex 



o 



In particular 

— X 

dx 6 = bF(e) where F(e) = 1 - e + C(e 2 ). (A.16) 
1 ~r 2ex 



Using (|A.16|) we conclude that (|A.13| ) and (|A.3|) (at n = 0) are identical on {c n } if 6 is 
chosen such that 

r 

6F(e) + = b(F(e) - 1) . (A.17) 

3=3 



We wish to find a solution to ( |A.17| ) that obeys ( |A.8| ), i.e. (from ( |A.16| )) one for which 



lim e ^o b = A. As A + Sj=3 1 = 0? such a solution exists, and takes the form 

6(e) = X(l + Ke + 0(e 2 )) (A.18) 

at small e where K is a number that may easily be determined. 

In summary, {c n } given by ( |A.15|) with b given by (|A.17|) , ( |A.18|) , solve ( |A.11|) for 
n > 1 and ( |A.3|) (with n = 0). {c n } therefore also approximately satisfy the true difference 
equations (|A.3|) for all n as long as |c n | 1 for all n > 1. But it is easy to verify that all 
\c n \ for all 7i > 1 are arbitrarily small at small enough e. Using the completeness of the 
Laguerre polynomials, 

I>»^ 2 ( (TTW< ft2 - (A ' 19) 



n=0 
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But 

"OO — x 



/■oo — x roo 

c =b dx >b dxe~ x (l-2ex) = 6(1 - 2e). (A.20) 

J 1 + 2ex J 

Combining flP9D and flOCD 

oo 

E c ' < 4e&2 ( A - 21 ) 

n=l 

establishing (|A.9| ) uniformly in n on our solution. Thus {c^} provides an approximate 
solution to the full nonlinear difference equations (|A.3|) for all n at small enough e. Fur- 
thermore, from ([A. 19 ) , this solution has finite energy. 



As {c n } obey the linearized recursion relation (|A.11[) and are small at small e, we 
can conclude, from the previous subsection, that d n takes the form (|A.7| ) for ne > 1. In 
order to estimate the behaviour of c n (e) for n C |- we formally expand the denominator 
in ( |A.15| ) in a power series in ex and integrate term by term, arriving at the asymptotic 
expansion 



oo 



X ^ / ^ \ 'til — I— Tl / n \ 'I I I 



|2 



( _ ir+ n (2er (A22) 

n!(m — nj! 



This expansion is useful only when the first few terms in the series in (|A.22|) are successfully 
smaller, i.e. for ne 1. 

A. 3. Generalization to Higher Dimensions 

In this subsection we will outline the generalization of the arguments of A.l and 
the construction of A. 2, for the case of the maximally isotropic noncommutativity in 21 
dimensions, i.e. a theory with noncommutativity matrix 0, all of whose eigenvalues are 
±i9. It is likely that these arguments can be further extended to generic O. 

We first note that a subset of the diagonal 9 = oo solutions ( [3.271) are (in non- 
dimensionalized coordinates) invariant under SO (21) rotations. These solutions take the 
form 



J 

Here 

*SV = EH a )= 2 '(- 1 )'4 l_1 V)- (A-24) 



where Lj ^(r 2 ) is an associated Laguerre polynomial. ( |A.24| ) is obtained from (|3.27| ) by 
repeated use of the identity 

n 

J2 L n-m^)Li(y) = L^+\x + y). 

m=0 
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D 



J + 1-1 



is a convenient normalization factor. 

When the noncommutativity matrix is maximally isotropic, the kinetic term in ( |2.3| ) 
is invariant under SO(2l) rotations of rescaled coordinates. Thus the corrections to an 
SO(2l) invariant 9 = oo solution, of the form (|A.23| ), are also SO(2l) invariant. 

Restricting to SO(2l) invariant functions, the arguments of section A.l are easily 
generalized. Any SO (21) invariant function takes the form 



n=0 



(A.25) 



The equation of motion implies that cj obey the following generalization of (|A.3| ) 



cj + J2 d 3 c 7 l = H( J + l - - (2J + l)cj + (J + l)c J+1 }. (A.26) 

J=3 

For large J (|A.26| ) and (|A.3|) are identical, hence all conclusions of section A.l carry over 
to this case. 

The perturbative construction of the solution that reduces to the SO (21) invariant 
Gaussian proceeds as in section A. 2 yielding the approximate result (good for small e) 



d, 



SDjT(l) Jo 



dx- 



x x l - 1 e- x L { j~ 1) (x) 



1 + 2tx 



(A.27) 
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